------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  ...\1903 JOP\trust_spd.log
  log type:  text
 opened on:  15 Mar 2019, 14:23:46

. 
. set seed X52dd25dd95a5bf6d0bef632d25b0dd210004410b

. 
. logit trust_spd pattern reform ///
>      saccounty sac_pattern sac_reform ///
>          black black_pattern black_reform

Iteration 0:   log likelihood = -1546.6054  
Iteration 1:   log likelihood = -1464.3242  
Iteration 2:   log likelihood = -1464.1688  
Iteration 3:   log likelihood = -1464.1688  

Logistic regression                               Number of obs   =       2328
                                                  LR chi2(8)      =     164.87
                                                  Prob > chi2     =     0.0000
Log likelihood = -1464.1688                       Pseudo R2       =     0.0533

-------------------------------------------------------------------------------
    trust_spd |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
--------------+----------------------------------------------------------------
      pattern |   -.486104   .1469679    -3.31   0.001    -.7741558   -.1980522
       reform |   .1850994   .1475611     1.25   0.210    -.1041152    .4743139
    saccounty |    .500561    .150845     3.32   0.001     .2049102    .7962118
  sac_pattern |   .3748702   .2151522     1.74   0.081    -.0468204    .7965607
   sac_reform |  -.0923874   .2202067    -0.42   0.675    -.5239846    .3392098
        black |  -1.457369    .260626    -5.59   0.000    -1.968187   -.9465514
black_pattern |    .166278   .3924985     0.42   0.672     -.603005    .9355609
 black_reform |  -.2051669    .375665    -0.55   0.585    -.9414567    .5311229
        _cons |   .4556613   .1020847     4.46   0.000     .2555789    .6557436
-------------------------------------------------------------------------------

. 
. estsimp logit trust_spd pattern reform ///
>      saccounty sac_pattern sac_reform ///
>          black black_pattern black_reform

Iteration 0:   log likelihood = -1546.6054
Iteration 1:   log likelihood = -1464.3242
Iteration 2:   log likelihood = -1464.1688
Iteration 3:   log likelihood = -1464.1688

Logistic regression                               Number of obs   =       2328
                                                  LR chi2(8)      =     164.87
                                                  Prob > chi2     =     0.0000
Log likelihood = -1464.1688                       Pseudo R2       =     0.0533

------------------------------------------------------------------------------
   trust_spd |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     pattern |   -.486104   .1469679    -3.31   0.001    -.7741558   -.1980522
      reform |   .1850994   .1475611     1.25   0.210    -.1041152    .4743139
   saccounty |    .500561    .150845     3.32   0.001     .2049102    .7962118
 sac_pattern |   .3748702   .2151522     1.74   0.081    -.0468204    .7965607
  sac_reform |  -.0923873   .2202067    -0.42   0.675    -.5239845    .3392099
       black |  -1.457369    .260626    -5.59   0.000    -1.968187   -.9465514
black_patt~n |    .166278   .3924985     0.42   0.672     -.603005    .9355609
black_reform |   -.205167    .375665    -0.55   0.585    -.9414568    .5311229
       _cons |   .4556613   .1020847     4.46   0.000     .2555789    .6557436
------------------------------------------------------------------------------

Simulating main parameters.  Please wait....
% of simulations completed: 11% 22% 33% 44% 55% 66% 77% 88% 100% 

Number of simulations  : 1000
Names of new variables : b1 b2 b3 b4 b5 b6 b7 b8 b9

. 
. 
. * Control group, non-Sac
. 
. setx median

. simqi, prval(1) genpr(control_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         0       median  
black_pattern |         0       median  
 black_reform |         0       median  
      pattern |         0       median  
       reform |         0       median  
  sac_pattern |         0       median  
   sac_reform |         0       median  
    saccounty |         0       median  
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~d=1) |   .6098958     .0245315      .561735    .6576745

Simqi generated the following new variable(s): control_pr

. 
. simqi, fd(prval(1) genpr(pattern_fd)) changex(pattern 0 1)

First Difference: pattern 0 1

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
         dPr(trust_~d = 1) |   -.118548     .0354168    -.1905204   -.0513339

Simqi generated the following new variable(s): pattern_fd

. simqi, fd(prval(1) genpr(reform_fd)) changex(reform 0 1)

First Difference: reform 0 1

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
         dPr(trust_~d = 1) |   .0451725     .0346325    -.0255881    .1109281

Simqi generated the following new variable(s): reform_fd

. 
. 
. simqi, fd(prval(1) genpr(sac_fd)) changex(saccounty 0 1)

First Difference: saccounty 0 1

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
         dPr(trust_~d = 1) |   .1112485     .0333656     .0410927    .1743506

Simqi generated the following new variable(s): sac_fd

. simqi, fd(prval(1) genpr(black_fd)) changex(black 0 1)

First Difference: black 0 1

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
         dPr(trust_~d = 1) |   -.339041     .0546986    -.4403484    -.233195

Simqi generated the following new variable(s): black_fd

. 
. 
. * Pattern group, non-Sac
. 
. setx pattern 1

. simqi, prval(1) genpr(pattern_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         0       median  
black_pattern |         0       median  
 black_reform |         0       median  
      pattern |         1         1     
       reform |         0       median  
  sac_pattern |         0       median  
   sac_reform |         0       median  
    saccounty |         0       median  
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~d=1) |   .4913479     .0252427     .4416748    .5400708

Simqi generated the following new variable(s): pattern_pr

. 
. 
. * Reform group, non-Sac
. 
. setx pattern 0

. setx reform 1

. simqi, prval(1) genpr(reform_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         0       median  
black_pattern |         0       median  
 black_reform |         0       median  
      pattern |         0         0     
       reform |         1         1     
  sac_pattern |         0       median  
   sac_reform |         0       median  
    saccounty |         0       median  
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~d=1) |   .6550683     .0248951     .6025325    .7030225

Simqi generated the following new variable(s): reform_pr

. 
. 
. * Difference in baseline probabilities
. 
. tabstat control_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
  control_pr |  .6098958  .0245315   .499419  .5705075  .5791652  .6101412  .6402872  .6498244  .7019231
--------------------------------------------------------------------------------------------------------

. tabstat pattern_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
  pattern_pr |  .4913479  .0252427  .4134671  .4504904  .4590231  .4906443  .5233883  .5329467  .5755939
--------------------------------------------------------------------------------------------------------

. tabstat reform_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
   reform_pr |  .6550683  .0248951   .573357  .6141338  .6230056  .6561877  .6854654  .6947679  .7314277
--------------------------------------------------------------------------------------------------------

. 
. gen dif_con_pat_pr = pattern_pr - control_pr
(1759 missing values generated)

. gen dif_con_ref_pr = reform_pr - control_pr
(1759 missing values generated)

. gen dif_pat_ref_pr = reform_pr - pattern_pr
(1759 missing values generated)

. 
. tabstat dif_con_pat_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_con_pa~r |  -.118548  .0354168 -.2362288 -.1790439 -.1661294 -.1176328 -.0749057 -.0610665 -.0180796
--------------------------------------------------------------------------------------------------------

. tabstat dif_con_ref_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_con_re~r |  .0451725  .0346325 -.0726962 -.0111401  .0020607     .0447  .0890543  .1005073  .1729228
--------------------------------------------------------------------------------------------------------

. tabstat dif_pat_ref_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_pat_re~r |  .1637205  .0363161  .0504028  .1023246  .1184993  .1646082   .210056  .2237759  .2704621
--------------------------------------------------------------------------------------------------------

. 
. 
. * Control group, Sac
. 
. setx median

. setx saccounty 1

. simqi, prval(1) genpr(sac_control_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         0       median  
black_pattern |         0       median  
 black_reform |         0       median  
      pattern |         0       median  
       reform |         0       median  
  sac_pattern |         0       median  
   sac_reform |         0       median  
    saccounty |         1         1     
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~d=1) |   .7211444     .0235628     .6725935    .7637218

Simqi generated the following new variable(s): sac_control_pr

. 
. simqi, fd(prval(1) genpr(sac_pattern_fd)) changex(pattern 0 1 sac_pattern 0 1)

First Difference: pattern 0 1 sac_pattern 0 1

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
         dPr(trust_~d = 1) |   -.022901     .0336147    -.0869007    .0421797

Simqi generated the following new variable(s): sac_pattern_fd

. simqi, fd(prval(1) genpr(sac_reform_fd)) changex(reform 0 1 sac_reform 0 1)

First Difference: reform 0 1 sac_reform 0 1

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
         dPr(trust_~d = 1) |   .0173248     .0345366    -.0528815    .0864405

Simqi generated the following new variable(s): sac_reform_fd

. 
. 
. * Pattern group, Sac
. 
. setx pattern 1

. setx sac_pattern 1

. simqi, prval(1) genpr(sac_pattern_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         0       median  
black_pattern |         0       median  
 black_reform |         0       median  
      pattern |         1         1     
       reform |         0       median  
  sac_pattern |         1         1     
   sac_reform |         0       median  
    saccounty |         1         1     
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~d=1) |   .6982433     .0248267     .6453392    .7448357

Simqi generated the following new variable(s): sac_pattern_pr

. 
. 
. * Reform group, Sac
. 
. setx pattern 0

. setx sac_pattern 0

. setx reform 1

. setx sac_reform 1

. simqi, prval(1) genpr(sac_reform_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         0       median  
black_pattern |         0       median  
 black_reform |         0       median  
      pattern |         0         0     
       reform |         1         1     
  sac_pattern |         0         0     
   sac_reform |         1         1     
    saccounty |         1         1     
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~d=1) |   .7384692      .025255     .6874039    .7865758

Simqi generated the following new variable(s): sac_reform_pr

. 
. 
. * Difference in baseline probabilities
. 
. tabstat sac_control_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
sac_contro~r |  .7211444  .0235628   .640891   .682765  .6899794  .7230167  .7505637  .7575682  .7965505
--------------------------------------------------------------------------------------------------------

. tabstat sac_pattern_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
sac_patter~r |  .6982433  .0248267  .6059438  .6572564  .6651274  .6994269  .7298627  .7384346  .7716174
--------------------------------------------------------------------------------------------------------

. tabstat sac_reform_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
sac_reform~r |  .7384692   .025255  .6663341  .6942319  .7041085  .7391102  .7705586  .7799336   .808377
--------------------------------------------------------------------------------------------------------

. 
. gen dif_sac_con_pat_pr = sac_pattern_pr - sac_control_pr
(1759 missing values generated)

. gen dif_sac_con_ref_pr = sac_reform_pr - sac_control_pr
(1759 missing values generated)

. gen dif_sac_pat_ref_pr = sac_reform_pr - sac_pattern_pr
(1759 missing values generated)

. 
. tabstat dif_sac_con_pat_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_sac~t_pr |  -.022901  .0336147 -.1255566 -.0794123 -.0665898 -.0234079   .020838   .032658  .0672836
--------------------------------------------------------------------------------------------------------

. tabstat dif_sac_con_ref_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_sac_co.. |  .0173248  .0345366 -.0887574 -.0410032 -.0259596   .015857  .0605629  .0719165  .1234329
--------------------------------------------------------------------------------------------------------

. tabstat dif_sac_pat_ref_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_sac_pa~r |  .0402259   .035973 -.0888581 -.0196232 -.0050002  .0405331  .0882016  .0993268  .1387143
--------------------------------------------------------------------------------------------------------

. 
. 
. * Control group, Black
. 
. setx median

. setx black 1

. simqi, prval(1) genpr(black_control_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         1         1     
black_pattern |         0       median  
 black_reform |         0       median  
      pattern |         0       median  
       reform |         0       median  
  sac_pattern |         0       median  
   sac_reform |         0       median  
    saccounty |         0       median  
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~d=1) |   .2708549     .0519685     .1763686    .3728518

Simqi generated the following new variable(s): black_control_pr

. 
. simqi, fd(prval(1) genpr(black_pattern_fd)) changex(pattern 0 1 black_pattern 0 1)

First Difference: pattern 0 1 black_pattern 0 1

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
         dPr(trust_~d = 1) |  -.0561169     .0708076    -.1918423    .0793682

Simqi generated the following new variable(s): black_pattern_fd

. simqi, fd(prval(1) genpr(black_reform_fd)) changex(reform 0 1 black_reform 0 1)

First Difference: reform 0 1 black_reform 0 1

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
         dPr(trust_~d = 1) |  -.0017425     .0757418    -.1426633    .1552992

Simqi generated the following new variable(s): black_reform_fd

. 
. 
. * Pattern group, Black
. 
. setx pattern 1

. setx black_pattern 1

. simqi, prval(1) genpr(black_pattern_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         1         1     
black_pattern |         1         1     
 black_reform |         0       median  
      pattern |         1         1     
       reform |         0       median  
  sac_pattern |         0       median  
   sac_reform |         0       median  
    saccounty |         0       median  
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~d=1) |    .214738     .0483838     .1319907    .3142481

Simqi generated the following new variable(s): black_pattern_pr

. 
. 
. * Reform group, Black
. 
. setx pattern 0

. setx black_pattern 0

. setx reform 1

. setx black_reform 1

. simqi, prval(1) genpr(black_reform_pr) listx

You have set the following values for the explanatory variables:

----------------------------------------
     Variable |    Value     Description
--------------+-------------------------
        black |         1         1     
black_pattern |         0         0     
 black_reform |         1         1     
      pattern |         0         0     
       reform |         1         1     
  sac_pattern |         0       median  
   sac_reform |         0       median  
    saccounty |         0       median  
----------------------------------------


Quantities of interest based on those explanatory values:

      Quantity of Interest |     Mean       Std. Err.    [95% Conf. Interval]
---------------------------+--------------------------------------------------
            Pr(trust_~d=1) |   .2691124     .0537131     .1817967    .3795863

Simqi generated the following new variable(s): black_reform_pr

. 
. 
. * Difference in baseline probabilities
. 
. tabstat black_control_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
black_cont~r |  .2708549  .0519685  .1109414  .1937225   .208166  .2680437  .3410606  .3603925  .4758728
--------------------------------------------------------------------------------------------------------

. tabstat black_pattern_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
black_patt~r |   .214738  .0483838   .086116  .1441427  .1572887  .2098373  .2795033  .2984121  .4463909
--------------------------------------------------------------------------------------------------------

. tabstat black_reform_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
black_refo~r |  .2691124  .0537131  .1310869  .1889597  .2029621  .2641011   .344798  .3643378  .4901089
--------------------------------------------------------------------------------------------------------

. 
. gen dif_black_con_pat_pr = black_pattern_pr - black_control_pr
(1759 missing values generated)

. gen dif_black_con_ref_pr = black_reform_pr - black_control_pr
(1759 missing values generated)

. gen dif_black_pat_ref_pr = black_reform_pr - black_pattern_pr
(1759 missing values generated)

. 
. tabstat dif_black_con_pat_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_bla~t_pr | -.0561169  .0708076 -.3072625 -.1705072 -.1443172 -.0555151  .0354018  .0559017  .1671732
--------------------------------------------------------------------------------------------------------

. tabstat dif_black_con_ref_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_black_.. | -.0017425  .0757418 -.2733926 -.1229407 -.0944477   -.00304  .0979903  .1295973  .2513745
--------------------------------------------------------------------------------------------------------

. tabstat dif_black_pat_ref_pr, s(mean sd min p5 p10 p50 p90 p95 max)

    variable |      mean        sd       min        p5       p10       p50       p90       p95       max
-------------+------------------------------------------------------------------------------------------
dif_black_p~ |  .0543744  .0698277 -.1796181 -.0528728 -.0298926  .0522122  .1439801   .172751  .3645485
--------------------------------------------------------------------------------------------------------

. 
. 
. drop b1-b9

. drop control_pr-dif_black_pat_ref_pr

. 
. log close
      name:  <unnamed>
       log:  ...\1903 JOP\trust_spd.log
  log type:  text
 closed on:  15 Mar 2019, 14:23:48
------------------------------------------------------------------------------------------------------------
